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We show that it is possible to entangle three different 
many-particle states by Bragg spectroscopy with nonclassical 
light in a Bose condensate of weakly interacting atòmic gases. 
Among these three states, two are of àtoms corresponding to 
two opposite momentum side-modes of the condensate; and 
the other is of single-mode photons of the output probe beam. 
We demonstrate strong dependence of the multiparticle en- 
tanglement on the quantum statistics of the probe light. We 
present detailed results on entanglement keeping in view of 
the possible experimental situation. 
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Entanglement is the most intriguing feature of quan- 
tum mechanics. Of late, it has been recognized as a 
wonderful resource for quantum information processing. 
Creation of various entangled states is the first step to- 
wards development of quantum communication. In re- 
cent times, many schemes for entenglement production 
have been proposed and demonstrated. Entanglement 
in a Bose condensate arises quite naturally. For in- 
stance, Bogoliubov's theory |ÏJ of Bose-Einstein coneden- 
sation (BEC) in weakly interacting gases predicts that, in 
the condensate ground state, a pair of partides moving 
with opposit momentum are absolutely correlated (en- 
tangled). This nonclassical feature makes Bose conden- 
sates of weakly interacting atòmic gases an excedent 
source for entanglement in massive partides. Many au- 
thors have proposed the production of nonclassical states 
|^-^| in a two-component Bose condensate. Bragg scat- 
tering has also been considered (ï(]] for producing entan- 
glement and squeezing. 

Entanglement between two partides is quite common, 
for example, EPR states, polarization states of twin- 
photons, down converted two-photon states in optical 
parametric oscillator and so on. In contrast, three par- 
ticle entanglement is not so common, though recently 
three-photon GHZ entangled states |11| have been exper- 
imentally realized. Because of two dominent momentum 
side-modes q and — q involved in the Bragg scattering, a 
condensate seems to be a suitable candidate for explor- 
ing tripartite entanglement among these two condensate 
side-modes and the probe photon mode. We demonstrate 
that Bragg scattering with two input light beams - one 
pump and the other probe, generates entanglement be- 
tween the output probe and the two side-modes excited 
due to the scattering of light. We analyze in detail the 
relevant parameter regimes where such tripartite entan- 



glement will show up. We consider three different input 
probes: coherent, vacuum and one-photon Fock state. 
Although vacuum field as an input has been considered 
earlier and the entanglement between the scattered field 
and one condensate side-mode (q) is discussed |ÏÏ|, here 
we report another interesting regime where strong entan- 
glement in long time limit persists. Our results show that 
although a vacuum or a coherent probe field can inducc 
entanglement between the scattered field and one excited 
side-mode (q) of the condensate, the other condensate 
side-mode (-q) remains immune to entanglement with 
the field. The most interesting result we obtain is that 
a nonclassical one-photon field state as an input probe 
can cause this other mode to become entangled with the 
field and thus generates the desired tripartite entangle- 
ment. We also find that nonclassical field is useful for 
detecting the entanglement between the two condensate 
side-modes. 

Thcre are mainly four types of processes in pump- 
probe or Bragg scattering of light by a condensate. First, 
one pump photon is transformed into a probe photon 
causing conversion of a zero-momentum atom into an 
atom of momentum q. In the second process, an atom 
moving with a momentum — q is scattered back into a 
zero-momentum atom by the transformation of a photon 
from pump to probe mode. The other two are the reverse 
processes of these two. Both the pump and the probe 
làser beams are detuned far off resonance from an elec- 
trònic excited state of the àtoms in order to avoid heating 
of the condensate. Because of interplay and simultainity 
of such processes, Bragg spectroscopy holds the key for 
generating many-particle entanglement in the motional 
states of the àtoms and photons. 

The Hamiltonian of the system H = Ha + Hp + Haf, 
with Hp = TïwiCj^Ckj + T1UJ2 Ck Ck 2 an d 
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where Ck(c k ) represents the annihilation(creation) oper- 
ator for a làser photon with momentum k, àk(à k ) is the 
annhilation(creation) operator for àtoms with momen- 
tum k and frequency u k = íl = (É 1 .d 13 )(É2.d 3 2)/ A 
is the two-photon Rabi frequency, where £1(2) i s the the 
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pump(probe) field amplitude with frequency uu^) an d 
the dij is the electrònic transition dipole moment between 
the states \i) and \j) of an atom. Here a s is the s-wave 
scattering length of the àtoms and V is the volume of the 
condensate. 

We assume that zero-momentum (k = 0) conden- 
sate state is macroscopically occupied and therefore the 
atom-atom interaction characterized by s-wave scatter- 
ing length for a weakly interacting atòmic gas is mainly 
due to the collision between zero- and and non-zero 
momentum àtoms. Applying Bogoliubov's prescription 
ào, à \A/Vo, whilc the condensate fraction Nq/N, with 
N being the total number of àtoms, and the number den- 
sity no = Nq/V remaining fixed in the thermodynamic 
limit, we convert the hamiltonian Ha into a quadratic 
form ]Ï2| . Further, we apply Bogoliubov's transformation 
jl2[ between particle and quasi-particle operators: àk = 
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is energy of Bogoliubov's quasi-particle. Here /i = £ m 
is the chemical potential with £ = (87moa s ) -1 / 2 be- 
ing a characteristic length scale known as the healing 
length. The condensate ground state energy is E g = 
\N n — ^J2k^o v t^k- We thus diagonalize Ha and 
rewrite the entire Hamiltonian in terms of Bogoliubov's 
quasi-particle operators àk- We consider the condensate 
ground state energy as the zero of the energy scale. By 
treating the pump light beam classically, the effective 
Hamiltonian can be written in the Standard form 
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where 5 = uj\ — u)2, q = ki — k2 and r\ = yN f q ÇÏ; where 
fq = u q — v q . In writing Eq.(4), we have retained only 
two dominent momentum side-modes of the condensate, 
and neglected all other modes under the Bragg resonance 
condition (S ~ u> q ). Equation (4) involves three cou- 
pled operators. The evolution operator ex.~p(—iH e fft/K) 
can not be disentangled into those of individual opera- 
tors, since there exists no such disentanglement theorem 
for such Hamiltonian. We therefore, emphasize that the 
Hamiltonian (4) can be solved exactly and more easily 
in Heisenberg picture. The Heisenberg equations of mo- 

tion for a triad of operators X = ^ct q àlq c£ 2 J can 

be written in a matrix form X = iuiq'M.X . The dynam- 
ics of these three coupled modes (two momentum side 
modes ±q plus one probe field mode) is controlled by 
the eigenvalues of M matrix. The long time behavior 
is oscillatory or growing depending on the momentum 



transfer q and the effective coupling constant 77. For 
a given value of q, if r\ exceeds a threshold value rjth, 
the long-time dynamics becomes hyperbolic dominated 
by the complex eigenvalues. In contrast, if r\ is less than 
77í/ l7 the long time behavior is oscillatory. Fig.l shows the 
dependence of rjth on the momentum q of the recoiling 
àtoms. In the phonon or quasi-particle regime (Çq < 1), 
rjth is much larger compared to that in particle regime 
(Çq » 1). For experimental situation of Ref. |ï^] with 
Na condensate, where quasi-particle regime (Çq ~ 0.47) 
is probed by Bragg scattering of light, r\th = 0.314kj^ 
with Lüq = (2-7r) x 4.7 kHz. The value of 77 used in that 
experiment is about 150 uj^ . 
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FIG. 1. The threshold value of r\ (rjth/Uq) as a function of 



To demonstrate many-particle entanglement among 
two motional (momentum) states of the condensate and 
one probe field state, following Ref. |Ï0| , |Ï4| ], we define 
two-mode entanglement parameter 

= ([A(rii - fij)] 2 ) I {{fii) + (fij)), i,j = q,-q,k 2 (5) 

where ((An) 2 ) = (tt, 2 ) — (n) 2 and the two particle 

number operators à^à q , à^ q à_ q and one photon num- 
ber operator c£ Ck 2 . The particle operators à are re- 
lated to the quasi-particle operators à by Bogoliubov's 
transformation. If £jj is below unity, the correspond- 
ing two- modes 'i', and 'j' are entangled. We specifically 
consider a condensate of 5 x 10 6 Na àtoms. First, we ex- 
amino entanglement by light scattering in phonon regime 
(£g < 1). Whilc choosing the relevant parameters, we fol- 
lowed the experiment of Ref. |Ï3[] , i.e., chemical potential 
u = 6.7 kHz, the pump and the probe fields intersecting 
at an angle 14° transfer a momentum q = 0.47£ _1 to the 
àtoms. Figure 2 illustrates the different two-mode en- 
tanglement parameters as a function of time with probe 
field being initially in a coherent state |/3) with average 
photon number \(3\ 2 — 1.0. Compared to long-time hy- 
perbolic limit (Fig.2(a)), the oscillatory limit (Fig.2(b)) 
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seems to be morè interesting from the two-mode or mul- 
timode entanglement viewpoint. As shown by Fig.2(b), 
entanglement between the (q, —q) and (q, fo) revives and 
disappears with time in the long-time limit, while in the 
hyperbolic limit, entanglement completely vanishes after 
a certain time (Fig.2(a)). Before vanishing, £ Çi fc 2 makes 
small oscillations eventually reaching a minimum (maxi- 
mum entanglement). As the coupling rj increases, these 
oscillations become less prominent (inset to Fig.2(a)), 
while for higher momentum (Çq > 1), these oscillations 
completely die down (inset to Fig.2(b)). The oscillations 
are caused by the interference of the three closely lying 
eigenvalues of M. From Fig.2, we infer that although 
the two condensate side-modes q and -q and one side- 
mode q and the scattered field mode k2 are entangled 
in different time and parameter zones, modes -q and k.2 
does not exhibit any mutual entanglement in any zone. 
However, by changing the character of the applied probe 
field, as we show below, it is possible to obtain mutual 
entanglement among all the three modes. By comparing 
Fig.2(b) with Fig.3(a) which exhibits £jj when the probe 
field is in vacuum, we further infer that Bragg scattering 
with either coherent or vacuum probe field generates al- 
most similar entanglement characteristics in the system. 
Although there arises mutual entanglement between the 
modes (q, — q) and (q, k/2) due to either of theses two 
fields, entanglement between (— q, k^) is not developed 
in either case. The results depicted in Fig.3(a) may be 
contrasted with those of Ref. [jï(J which studied entan- 
glement between q and k.2 modes in the hyperbolic limit 
for a vacuum probe. In contrast, the oscillatory limit we 
describe in the present paper is a new result not discussed 
before thus far. We emphasize that entanglement in the 
oscillatory limit is more significant; since it sustains at 
long time, while in the hyperbolic limit it does not exist 
at long time. We also point out that these two limits 
depend on the two parameters q and ry, and not due the 
dynamical processes of the three coupled modes. 

Next, we consider the probe field as a one-photon 
Fock state |^]. In this case, all the three modes exhibit 
mutual entanglement at different times and duration as 
is evident from the Fig.3(b). Thus three-mode or tripar- 
tite entanglement can be generated in a Bose conden- 
sate by Bragg spectroscopy with nonclassical light fields 
(Fock states) which constitutes another central result of 
the present paper. It should be pointed out that neither 
vacuum nor coherent probes are able to generate such 
tripartite entanglement in any parameter space. Why 
one-photon field state and neither vacuum nor the coher- 
ent state can generate tripartite entanglement may be 
explained by examining the respective phae-space den- 
sity distribution functions. Since, coherent state is a dis- 
placed vacuum state, both of them have almost similar 
phase-space structure, except a displacement of the equi- 
librium position. We have developed a formalism based 
on time-dependent Wigner distribution functions which 



provides significant insight into the entanglement prop- 
erties in various cases discussed in this paper. 
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FIG. 2. (a) Two-mode entanglement parameter Ç 9 ,- 9 (sòlid 
line), Ç q ,k 2 (dashed lines) and Ç- 9 ,fc 2 (dotted line) as a func- 
tion of time t in microsecond (ps). Initially, the condensate 
is assumed to be in the ground state and the probe field is 
in a coherent state with average photon number (n p ) = 1. 
The other parameters: q = 0.47£~\ uig = (2tt) x 4.7 kHz, 
the two-photon Rabi frequeney fi = 8Hz, rj/íjq = 0.34 cor- 
responding to hyperbolic regime, since r}th/Wq = 0.313. The 
eigenvalues of M matrix are 1.07975, - 0.69788 + 0.14153Í 
and -0.69788 - 0.14153Í. (b) Same as in (a), but fi = 7 
Hz corresponding to oscillatory regime (rï/uil 3 = 0.29) with 
eigenvalues 1.06235, -0.79248, -0.58587. The inset to Fig.(a) 
shows the variation of Ç 9 ,fc 2 with t in /_ts for fi = 16 Hz and 
q = 0.47Ç" 1 . The inset to (b): £, qM Vs t in fis for fi = 16 Hz, 
q = 2C 1 - 

Now, to know the entanglement property in the par- 
ticle regime (£9 >> 1), we show the £y as a function 
of time for q = 8.329£ _1 in Fig.4. While considering 
phonon (oj^ oc q) or particle ( oj^ oc q 2 ) regime, exact 
value of ujq should be taken, otherwise an approximation 
in Lúq may lead to wrong results if 77 is close to the thresh- 
old value r]th- Here also we consider a Na condensate of 
5 x 10 6 àtoms, but the other parameters (/i, q) are chosen 
from the experimental paper of Ref. JÏ6| From Fig.4(a), 
we observe that before reaching the hyperbolic regime, 
all the three modes remain entengled for a considerable 
duration. The inset to Fig.4(a) shows the variation of £y 
as a function of the two-photon Rabi frequeney íl at time 
t = 10/ís. From this figure we conclude that the tripar- 
tite entanglement is signicant at weak coupling (low Í2 
or low rf) near the threshold (77^). Figure 4(b) exhibits 
entanglement property in the same particle regime but 
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in the long-time oscillatory limit. In this case, unlike the 
other cases, entanglement among the three modes per- 
sists for almost all the time, although in the long-time 
limit the modes q and — q (sòlid curve) are marginally 
entanglcd. But, in contrast, the modes q and k2 remain 
largly entanglcd for all the time. 
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FIG. 3. (a) Samc as in Fig.2(b), but for the probe field 
being initially in vacuum (cavity). The inset to Fig.3 (a) 
shows the Q parameter Q p of the output probe field as a 
function of time in /is. This Q parameter always remains 
greater than unity. (b) Same as in (a), but for probe field 
being initially in one photon Fock state. The inset to Fig.3 
(b) shows the corresponding Q p as a function of time in fj,s. In 
this case, the Q parameter becomes less than unity implying 
nonclassicality of the output light field. 

To see whether the photon number fluctuation of the 
output (scattered) light field has any connection with 
the entanglement parameters discussed so far, we calcu- 
late the Mandel Q parameter Q p = ((Ank 2 ) 2 )/(nk 2 ) We 
find that neither a coherent nor a vacuum field as an in- 
put probe can generate nonclassical light output, while a 
nonclassical one-photon Fock state can do so. The time- 
evolution of the Q parameter has a strong link with the 
time-evolution of the entanglement parameter i.e., 
entanglement between the atòmic modes q and — q. With 
a nonclassical input probe, both the parameters Q p and 
£q,-q become less than unity at the same time (Fig.3 (b) 
inset) indicating that the entanglement between the two 
atòmic side-modes may be inferred by measuring Q p us- 
ing beam splitter at the output of the probe beam and 
detecting the photon number fluctuations with two de- 
tectors. 

In conclusion, we have presented an exact treatment 



of the dynamics of excitations in a condensate by Bragg 
scattering of light under Bragg resonance condition. We 
have shown that by using a nonclassical light (Fock 
states) as input probe, it is possible to generate tripartite 
entanglement in the many-particle states of three mo- 
tional modes which inelude two momentum side-modes 
of the condensate and one probe field mode. The rel- 
evant physical parameter regimes where such tripartite 
entanglement can be observed are analyzed. The im- 
portant question which needs detailed study is how to 
detect the generated entanglement. Our results suggest 
that by measuring the second order correlation function 
of the scattered field, it is possible to measure the en- 
tanglement between the two modes of the condensate. 
The number variance of the q and q momentum àtoms 
may also be measured by outcoupling them with a change 
in the hypcrfine spin state by applying Doppler-sensitive 
two radio frequeney pulses and then by measuring the 
intensity variation at the output of the two pulses. 



0.001 



1 

0.01 
0.0001 




1 10 100 
£2(Hz) 




„ Vv^W^' (a) 




aF 0.01 



0.0001 



10 

t(jis) 

FIG. 4. (a) Same as in Fig.3(b), but (j, = (2tt) x 1.23 kHz, 
q = 8.329£ -1 ,a;f = (2tt) x86.65fcHz, fi = 7 Hz, jj = 0.0285wf 
and rj th = 0.0071wf . The eigenvalues are 1.00041, -0.99315 + 
0.0277H, -0.99315 - 0.0277H. (b) Same as in (a), but fi = 1 
Hz corresponding to oscillatory regime (r/ = 0.0041wf ) with 
eigenvalues 1.00001, -0.9987, -0.9872. The inset to Fig.(a) 
shows the variation of £jj with two-photon Rabi frequeney fi 
in Hz at time t = 10/ís with other parameters remaining the 
same as in the main figure of (a) . 
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